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CALCULATION OF UNSTEADY TRANSONIC AERODYNAMICS 
FOR OSCILLATING WINGS WITH THICKNESS* 

By S . Y . Ruo and J . G . Theisen 
Lockheed -Georgia Company 


SUMMARY 


An analytical approach is presented for determining some of the nonlinear characteristics 
of the transonic potential flow equation for finite thickness wings undergoing oscillatory 
motion at sonic flight speed in an inviscid, shock-free fluid. The thickness effect is 
accounted for in the analysis through the use of the steady local Mach number distribution 
on the wing, and corresponding mean flow parameters over the wing are considered to be 
constants in a small finite region. This latter assumption, equivalent to the concept of local 
linearization, permits the non linear differential equation for the velocity potential to be 
reduced to a linear one with variable coefficients containing the local Mach number. By 
means of an appropriate coordinate transformation, the equation becomes identical with the 
usual linearized transonic unsteady-flow equation with constant coefficients. Numerical 
results for wings with finite thickness can then be obtained by any computational technique 
that is suitable for use with the transonic linearized theory. For demonstration purposes 
herein, the well-known sonic -box computer program is adopted. 

Stability derivatives and generalized aerodynamic force coefficients are shown for delta, 
rectangular and parabolic wing planforms over a range of reduced frequencies, with wing 
thickness and aspect ratio as parameters. Comparisons are made with available linearized 
theories and test data in the transonic flow regime, and reasonable agreement in trends is 
indicated. Relatively small thicknesses of 5 percent can have significant effects on gener- 
alized aerodynamic forces, e.g., on the order of 25 percent for rectangular wings of low 
aspect ratio, although the results for delta wings are much less. Still greater thickness 
levels may show a reversal of magnitude trend, according to these results, with corresponding 
effects possible for flutter instabilities. 


* Portions of this report were included in AIAA Paper 73-316 which was presented at the AIAA 
Dynamics Specialists Conference, Williamsburg, Va., March 19-20, 1973. 



INTRODUCTION 


'The study of flutter and other aeroelastic responses of an aircraft requires adequate 
knowledge of the forces acting on three-dimensional wings in oscillatory motion. Such 
aeroelastic problems are frequently critical in the transonic speed range. The physical prob- 
lem is governed by a nonlinear partial differential equation with nonlinear boundary condi- 
tions, for which an exact solution is not known to exist. The basic small perturbation 
equation governing the velocity potential for transonic flow over a thin wing or a slender 
body is well known (e.g., refs. 1 and 2). However, the nonlinearity which remains present 
prevents closed-form solutions from being obtained except in a few special cases (ref. 1). 

For a low-amplitude, high-frequency oscillation where the unsteady part is considered to be 
a small disturbance to the steady part, the steady-state properties can be completely 
uncoupled from the unsteady equation, and the governing equation for unsteady transonic 
flow can be linearized (refs . 2 and 3) . 

Almost all unsteady transonic flow theoretical work lies within this framework of 
linearized theory where the thickness effect of the wing is neglected. However, oscillatory 
transonic aeroelastic instabi lities of lifting surfaces (i.e., flutter) often occur at frequencies 
below the range of validity of the transonic linearized theory, thus indicating a need to 
account for at least the more prominent effects of finite thickness. Recent studies (ref. 4) 
show that wing thickness, which enters into the mathematical nonlinearity, can significantly 
affect unsteady pressures in the supersonic range also. The object of this work, therefore, is 
to develop an appropriate method of accounting for the more important effects of finite wing 
thickness in order to predict transonic oscillating aerodynamic parameters at frequencies 
lower than those for which linearization is valid. The present study is limited to attached, 
shock-free flow . 

An important consequence of the previously mentioned linearization is the suppression 
of deviations in local Mach number from freestream value. Since these deviations have 
appreciable effect on propagation of pressure disturbances over the lifting surface, signifi- 
cant improvement in the theory may be accomplished by "recoupling" the steady and 
unsteady flow parameters in a manner that at least accounts approximately for variations in 
mean local Mach number caused by finite wing thickness. In the present development, this 
is achieved by considering all of the steady-flow parameters over the wing to be invariant 
within a small finite region . This latter assumption, equivalent to the concept of local 
linearization, permits the nonlinear differential equation for the velocity potential to be 
reduced to a linear equation with variable coefficients containing the local Mach number. 

By means of an appropriate coordinate transformation, the equation becomes identical with 
the usual linearized transonic unsteady-flow equation with constant coefficients. Hence, 
numerical results for a lifting surface in the transformed space may be obtained by any com- 
putational technique that is suitable for use with the transonic linearized theory. For 
demonstration purposes herein, the well-known sonic-box computer program (refs. 5 and 44) 
is adopted . 
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Since this investigation is exploratory, the limitations on wing planform geometry im- 
posed by the sonic-box method of reference 5 have not been removed. The types of wings 
which can be treated are those having unswept trailing edges without control surfaces. Also, 
because of the transformation method used, the steady flow everywhere over the wing must 
not be very different from that of the freestream. The mean angle of attack, twist and cam- 
ber of the wing are considered to be small, and the Mach number distribution on the wing is 
assumed to depend only on the wing thickness. Therefore, the mean local Mach number on 
a corresponding point of the upper and lower wing surfaces is considered to be the same in 
the present work. However, the extension of the present method to wings with different Mach 
number distributions on the upper and lower surfaces is self-evident; 

Calculations using the present method were made to evaluate the thickness contribution 
to the unsteady aerodynamic forces by comparison with cases not having any thickness effect. 
The wings considered were: (a) rectangular wings with a bicircular-arc profile; (b)delta 
wings with an elliptic lateral cross-section; and (c)parabolic wings with a double-wedge 
lateral cross-section. Comparisons are made with the sonic box method, for cases without 
thickness effects, as well as with available test data on finite thickness wings. 


A detailed description of the computer program is available as a supplementary report 
(refo 6). 


SYMBOLS 


A 

a 


B 

b 

c 

C L 

a 

c M£ +C M q 

C P 

D/Df 

E 

e, exp 

f 

H 


aspect ratio or influence coefficient 

dimensionless pitch axis location measured positive aft from most forward 
point on wing (reference length = b) 

B(x,y,z,t) = 0 defines the dimensionless body surface at any instant 
(reference length = b) 

reference length - wing centerline chord (dimension = L) 

local speed of sound (dimension = L/T) 

lift-curve-slope 

damping-in-pitch 

pressure coefficient 

substantial derivative 

box approximation of the wing planform 

exponential function 

dimensionless mode shape 

dimensionless length of box side (reference length = b) 
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dimensionless distance from z — 0 (reference length - b) 


'/ 1/ k 

Im 

k 

L 

L.. 

■I 

M 

& 

q 

Re 

S 

s 

T 

t 

U 

00 

U, V,W 

W 

*,y,z 

a 

r 

Y 

V 


unit vectors in x-, y-, z-direction, respectively 
imaginary part of the complex quantity 
reduced frequency, uub/U^ 
unit of length 

generalized aerodynamic force coefficient 
local Mach number 
reads "order of" 


/ 2 2 2 

local flow speed, \ u + v + w (dimension = L/t) 

real part of the complex quantity 

dimensionless planform area for full wing (reference area = b^) 
dimensionless local semi -span (reference length = b) 
unit of time 

dimensionless time (reference time = b/U^ 
freestream velocity (dimension = L/T) 

dimensionless velocity components in x-, y-, z-direction, respectively 
(reference velocity = U^) 

wake 

dimensionless cartesian coordinates (reference length = b) 
angle of attack 
acceleration parameter 
ratio of specific heats 

i ^- + i + / dimensionless gradient operator 
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2 

V 

6 

0 .. 

>1 


a 


$ 

<P 

cp 

a 


oj 


— *-4- — *-4- — ^ , dimensionless Laplacian operator 
c* dy^ Sz^ 

dimensionless amplitude, see equation (60) 

phase angle of the generalized aerodynamic force coefficients, (dimension 
degree or radian) 

maximum semi-span-to-chord ratio 
thickness-to -chord ratio 

dimensionless small perturbation velocity potential, $ = 0+ cp 
steady part of $ 
unsteady part of $ 

— 2 

total velocity potentia I, 0 - U^x + 1) (dimension - L /T) 
angular velocity (dimension = radian/T) 


Subscripts 


1,2 


d 


1 r I 


*/y,z 

xy,xt 


denotes the magnitude of the oscillatory quantity 
denote steady and unsteady conditions, respectively 
denotes conditions of the freestream 
denotes doublet solution 

denote summation indices in x and y directions, respectively 
denote first-order partial derivatives 
denote second-order partial derivatives 


Superscripts 

4- denotes upper surface 


denotes lower surface 



denotes sonic condition 


* 


denotes dimensional quantity or operator 

(i.e., ( ),. . , = ( ),. . | x (reference quantity)) 

v v dimensional 'dimension less n 

denotes vector quantity 

denotes transformed quantity 

denote first and second total differentiations 


UNSTEADY TRANSONIC FLOW EQUATIONS AND SOLUTIONS 


The Basic Equation for Potential Flow 


The exact, nonlinear, unsteady flow equation to be satisfied by the velocity potential 
Q , as given in reference 7, is 

V 2 fi - -L (fi) + (q 2 ) + q • v IX q 2 ) 

c L&t at 

where q = Vfl and c is the local speed of sound. All quantities in equation (1) are dimen- 
sional . 



Small Perturbation Potential Equation 


In order to obtain the small perturbation velocity potential equation, the small pertur- 
bation potential $ is introduced by writing 


^(x/yfZft) - U^Lx + $(x,y,z,t)], 
then, the velocity components are expressed as 

u = n-= u (i +I-) 1 

X °° v x 

v = fi-= u ¥- ► 

y 00 y 

w = n-= u §- 

z “ z 


( 2 ) 


(3) 
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The substitution of equation (3) into equation (2), with the introduction of the following 
set of dimensionless variables 


U t 


$ 


* = E ' y = l ■ Z = E ' t = TT ' and 1 "E 


results in the expression. 


2 - U » 


U 


v $ — *- [$ + 2$ + $ ] = {2($ $ +$$ +$$) 

2 xx xt tt 2 ' x xf y yt z zt' 


+ 2 ($ $ +$$ +$$ ) 

x xx y yx z zx 

+ $($$ +$$ +$$ ) 
x x xx y yx z zx 

+ §($<$ +$$ + $ $ ) 
y y xy y yy z zy 


+ $($$ + $ $ + $ $ )} # 

z x xz y yz z zz 


(4) 


Common ly, the perturbation is assumed to be so small that the quadratic and cubic terms on 
the right-hand side of equation (4) may be neglected in comparison with the linear term on 
the left-hand side. Thus, equation (4) simplifies to 


U 


U* 


U +$ +$ 7T (2$ + $ ) = 0 


xx yy zz 


xt tr 


(5) 


This may be considered as a linearized equation with variable coefficients. A further sim- 
plification is usually made by replacing the local speed of sound with that of the free- 
stream by employing the Bernoulli equation: 


9 2 

2 c 

c 03 


1 


—5 5 - (Y - 1) k+* +iW 

U U L t x z 

00 00 

Thus, equation (5) becomes the well-known differential equation of acoustic theory 

(1 $ +§ +$ - 2M 2 $ -M 2 $ =0. 

** xx yy zz 00 xt 00 tt 


( 6 ) 


(7) 


The equation of motion has been considerably simplified by using the small perturbation 
method; but this linearized equation does not always yield appropriate solutions for some 
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physical flow phenomena in transonic flow. In the study of transonic flow a nonlinear term 
is normally retained to account for the nonlinear effects which usually occur in this flow 
regime. Landahl (ref. 2) derived the small perturbation potential equation for unsteady 
transonic flow by retaining the first- and second-order term contribution from the 
x-derivative term in the continuity equation and only the first-order term contribution from 
the y-, z-, and t-derivative terms. The resulting equation is 

(i - m3 -m3<y+ o« x 

If the non-constant coefficient in equation (5) is approximated by 


$ + $ + $ - 2M 2 § - M 2 § = 0 . (8) 

xx yy zz ® xr co tt 



2 u 9 

= H_^M 2 (l-2i x ) 

c q 


and only the first order terms are retained, equation (5) becomes 

(1 -M 2 )? +§ + $ -M 2 (2§ + * ) = 0 . (9) 

' xx yy zz ' xt tt 

This equation maintains the mixed flow property usually observed in the transonic flow 
regime. It is not, however, in the form commonly used in transonic flow studies because of 
the presence of non -constant coefficients. Use of the freestream Mach number may be 
physically justified in place of the local Mach number in equation (9) as long as the local 
Mach number is not very different from unity. That is, both equations (7) and (9) are 
acceptable based on the analysis, but equation (7) does not preserve the mixed flow 
properties of the transonic flow problem. Thus, instead of equation (7), equation (8) is 
being used more frequently. 

Both equations (8) and (9) are capable of representing mixed flow but neither equation 
is easily tractable as required to satisfy the boundary conditions. These include the require- 
ments that the small perturbation velocity vanishes at infinity and the flow must be tangent 
to the solid surface. Also, in unsteady flow studies, the equation is further simplified and 
linearized by assuming that the unsteady potential is a small perturbation over that of the 
steady part. 


Unsteady Small Perturbation Potential Equation 


Under the assumption that the unsteady perturbation potential is small relative to that of 
the steady part, the steady and unsteady parts in the small perturbation equation may be 
completely separated when the unsteadiness (e.g., the frequency of oscillation) is high. 
Subsequently, the unsteady form of the equation is usually linearized. 
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Constant Coefficient Equation. - Following Landahl (ref. 2) the small perturbation 
potential $ may be written as 

§(x,y,z,t) =0(x / y,z) + cp(x,y,z,t) . (10) 

Then, for the case of oscillatory motion where the amplitude is much smaller than the thick- 
ness ratio, equation (8) becomes 

(1 - M 2)d> + 0 + 0 = M^(y + 1)0 0 (1 1) 

v oc/t'xx *yy r zz 20 ~x ^xx 7 


and 


0 - m3 


2 

CD + CD + CD - 2M CD 

xx yy zz 00 xt 


M 




( 12 ) 


where equations (11) and (12) correspond, respectively, to steady and unsteady transonic 
flows. Equation (11) involves only the steady flow quantities but equation (12) contains a 
coupling term between the steady and unsteady quantities. The term cpxCp xx bas been 
ignored because cp is considered to be a perturbation on 0 , so cpxCp xx is a higher order 
term than cp 0 or cp 0 . 

X r XX XX r x 

Landahl (ref. 2), following the fundamental development of Lin, etal. (ref. 8), showed 
that with the assumptions that cpj. = Gf (kcp x ) and k » 1 0 X I everywhere in the flow, the 
coupling term on the right-hand side of equation (12) can be neglected. Thus equation (12) 
becomes linearized and completely decoupled from the steady flow quantity, as follows: 


(1 - M + cp -i- cp - 2M^ cp — cp = 0 . (13) 

^ xx yy zz 00 xt °° tt 

Since in transonic flow 1 - ( 0^) , equation (1 3) may be further simplified to 

2 2 

cp + cp - 2N\ cp - M cp = 0 . (14) 

yy zz 00 xt °° tt 1 


Equation (14) invalid when 


i -m| sef(^) . 


In reference 9, Miles presented an excellent discussion of the conditions for the 
simplification of the small perturbation equation in various flow regimes. 

The freestream Mach number M^ in equation (14) can be eliminated by a transforma- 
tion using Y = M ro y and Z = MaoZ so that only the case when Moo = 1 needs to be considered. 
That is, for sonic flow, 
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( 15 ) 


cp +cp - 2p - cp = 0 . 

yy zz xt tt 

When the wing is performing oscillatory motion, the unsteady small perturbation 
potential can be expressed as 


, ,\ / \ i co t , % i kt 

c P( x /y/ z /l') =c P 0 ( x ,y,z)e =cp o (x,y,z)e 


( 16 ) 


and equation (14) becomes 


2 

cp + cp - M k (2icp - kcp ) = 0 . 

o “ v o o 

yy zz x 


(17) 


Non-constant Coefficient Equation.- If the steady local Mach number M is taken to be 
the mean local Mach number, equation (9), with the relation in equation (10), can be 
easily written as 


(1 - M 2 )^ 


+ (b + 

xx yy 


<t> = o 

r zz 


(18) 


and 


(1 - M^)cp + cp^ + cp^ - M 2 (2cp^ + cp t4 ) = 0 , 


XX 


zz 


r xt r tt' 

where equations (18) and (19) govern, respectively, the steady and unsteady flow. 


(19) 


By following the same reasoning used in obtaining equation (14) a similar equation can 
be derived for unsteady flow with non-constant coefficients, namely 


cp + cp 

yy 


- M 2 (2<P xt + cp n ) = 0/ 


( 20 ) 


or 


? 

cp + cp - M k(2i cp - k cp ) = 0 . (21) 

o o o o 

yy zz x 

Equation (21) is the basic equation considered in the present study, to investigate thickness 
effects on the aerodynamic forces, in which the steady local Mach number is treated 
as a parameter. 
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Boundary Conditions 


Tangency Condition.- For an inviscid fluid flowing over a solid surface, it is assumed 
that at the surface the normal component of the relative velocity is zero and slip may occur; 
there is no boundary layer separation and no shock formation. Thus the tangency condition 
can be written as (ref. 10) 


D B 
D T 


= 0 , 


( 22 ) 


where B(x, y, z, t ) = 0 defines the body surface position at any instant, and 

D S 

— = — + u + V — + w 

Dt Bt Bx By Bz 

For a small perturbation flow, equation (22) can be written in dimension less form as 

B + (1 + « )B + « B + « B =0, on B = 0 . (23) 

t x x y y z z 

To the first order, for a general thin body (§ « 1), equation (23) can be approximated by 

B + B + $B + $> B = 0, on B = 0 , (24) 

t x y y z z ' x ' 

and for a planar body, equation (24) can be further simplified if only the first order terms 
are retained: 


B + B + $ — 0, on z — 0 . (25) 

t x z x 7 

Equation (25) is the linearized tangency condition for the planar case. 

Fora thin body performing a small oscillation from its mean position, the instantaneous 
position of its upper surface may be approximated by (fig. 1) 

z = ^(x/y) + h 2 (x,y,t) 

or 

B(x,y,z,t) = z - hj(x,y) - h 2 (x,y,t) . (26) 

By substituting equation (26) into equation (24), one obtains 
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( 27 ) 


$ = h. + (h 0 + h« ) + $ (h, + h 0 ) , on B — 0 . 

z 1 2 2 l y I 2 

x x t ' y y 

Since the body under consideration is thin and its departure from the mean position is small, 
the spatial derivatives of hj and h 2 may be considered to have the same order of magni- 
tude and $y is much less than unity; thus, equation (27) may be linearized (i .e. , eq. (25)). 
When equation (10) is introduced in the left-hand side of the linearized form of equation 
(27), one arrives at 


and 


*z = h l 


cp z =h 2 + h 2 


> on z = 0 , 


(28) 


where the steady and unsteady tangency conditions have been uncoupled. 

For a body performing oscillatory motion 

^(x^y/t) = h 0 ( x /y) elllJ ' r= h Q (x,y)e lkf , 

the linearized unsteady tangency condition in equation (28) may be written as 


cp = h + ikh , on z = 0 . 
o o o' 

z x 


(29) 


Pressure Coefficient .- The linearized pressure coefficient written in terms of the small 
perturbation potential is 


C p (x,y,0,t) = - 2($ f + $ x ). (30) 


Substitution of the relationship in equation (10) into equation (30), with C p written as a 
linear combination of the steady and the unsteady parts, results in 


C p ( x /YfO) “ 

(steady part) 

(31) 

and 



C p 2 ^ X,y ' 0/t ^ = ~ ■ 2 ( CP x + Cp P 

(unsteady part) . 

(32) 

If equation (16) is substituted in equation (32) the following expression is obtained 
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(33) 


where 


C = - 2(cp + i k cp ) 

p o o ' 

■ rs v 


c p ( x ,y,0) =c p (*/y/0/t) e lkt . 


For locations in the wake of a wing, = 0, or 

cp + i k cp = 0 . 
o o 

x 


(34) 


Forms of the Nonlinear Term in Steady Flow 


In order to maintain the mixed-type flow in the study of steady transonic Flow, the 
variable coefficient of the first term in equation (18) is often approximated by the general 
form 


(1 - M 2 ) -0 - m3 -K 




(35) 


where K, which represents the coefficient of the nonlinear term, assumes several forms as 
listed below: 


= (Y+ 1)M l 

(e.g . , 

ref. 11 or eq . (1 1)) 

(36) 

= (Y+ 1) 

(e.g.. 

refs. 12, 13, 14, 15) 

(37) 

= M 2 [2+(y- 1)M 2 ] 

(e.g.. 

ref . 1 6) 

(38) 

2 

1 -M 

CO 

(e.g ., 

refs. 17, 18) 

(39) 


+ 


V + 1 ’ MJY+ 1) 

K (Y + 1 ) 0 (X*) 

d> 00 r xx 

~ XX 


K (Y+ 1) [B"(x*) • 4n t] ( X -X*) 


(e.g . , refs. 19, 15, 20) 

(e.g . , ref. 21) . 

.2 


(40) 

(41) 


Equation (38) can be obtained by expressing M in terms of and utilizing the 
Bernoulli equation relating the local and freestream speed of sound, and retaining the first 
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order terms. When Moo tends to unity, equations (36), (38), and (39) all reduce to equation 
(37) . Equation (39) was obtained from the mass-flow consideration . Equations (40) and 
(41) were obtained more or less intuitively. Equation (40) forms the basis for the parabolic 
method. Equation (41) was needed to treat the body of revolution for the case of M® = 1 
only; the form of K would be different when the freestream Mach number is not unity. 

When the wing is of very small aspect ratio and 1 - M = & (0 ), K may be set to 
zero (e.g., refs. 2, 22). 


Theoretical Methods to Obtain Numerical Results 


In all the practical unsteady transonic flow studies to date, only the sonic freestream 
flow case has been considered. This is because the linearized equation with constant coef- 
ficients, namely, equation (13) has been used and a simple transformation can be used to 
eliminate M® in the equation. Then a similarity rule can be applied to obtain information 
for those cases where M» is slightly different from unity. The linearized theory has been 
used with some success in comparison with test data. But due to the possible presence of 
wind tunnel wall interference, shock formation, and flow separation it is sometimes difficult 
to make a direct comparison of the test data with the theoretical results and to identify the 
source of disagreement. Because of this, only the theoretical methods leading to numerical 
results are mentioned in the following paragraphs and no attempt is made to evaluate the 
validity of each method. 


Linearized Method.- Mangier (ref. 23) used the pressure doublet distribution to treat 
the delta wing at sonic speed, for small k, by letting M a approach unity from the subsonic 
solution. Miles analyzed the low aspect-ratio rectangular wing at sonic speed in reference 
24. Later, Runyan and Woolston (ref. 25) solved the case of an oscillating rectangular wing 
at sonic speed by extending the kernel function method for general planforms in subsonic 
flow. In reference 2, Landahl extensively treated the delta and rectangular wings by 
solving the boundary value problem in a Fourier transformed plane. Using the delta wing 
solution as a tool, Landahl also studied low aspect-ratio wings, wing-body combinations, as 
well as wings of general form composed of straight segments. Davies (ref. 3) used lifting 
surface theory to treat general planforms for large k. An elliptic cone performing plunging 
oscillations at sonic speed was investigated by Albano and Andrew (ref. 26) for the effects 
of blunt leading edges and wing thickness. 

Rodemich and Andrew (ref. 5) developed a sonic -box method for non-swept trailing 
edge wings with arbitrary wing planforms; and, later, Stenton and Andrew (ref, 27) 
extended the method to include swept trailing edges with control surfaces. 

Parabolic Method .- In the linearized method the unsteady flow part has been com- 
pletely uncoupled from the steady part; that is, the thickness effect has been neglected. An 
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attempt- to take the coupling of the stead/ part Into account, at least to the first order, has 
been made by Telpel (ref. 28) for two-dimensional flow and Liu et al. (refs. 29, 30) for 
axisymmetrical flow at sonic speed. This is based on the parabolic method concept (eq. 
(40)) introduced by Oswatitsch and Keune (ref. 19) for steady axisymmetrical flow. An 
equation in the form of equation (12) is used with the term cp cp ignored and with cp cp 
written as rcp^/[/V\£(y + 1)]. 


Proposed Method.- It has been shown by Kacprzynski, et al. (ref. 4), for supersonic 
flow that thickness^effects on the lift, aerodynamic center, and flutter speed of a delta 
wing could be very significant. It is not difficult to conceive that the thickness could also 
have a marked effect on the aerodynamic forces on a finite wing over the transonic flow 
range. The thickness effect is introduced through the use of the local Mach number of the 
wing at its mean position in the unsteady, small perturbation equation. 

The proposed equation to be used in this investigation is equation (20) in which M is 
the steady local Mach number representing the effect due to wing thickness. The solution to 
equation (20) satisfying the required boundary conditions is obtained by using the concept of 
local linearization attributed to Spreiter and Alksne (ref. 31) as described in the following 
section . 


LOCAL LINEARIZATION 


Local Linearization Concept 


The thickness of the airfoil, the formation of shocks, the interaction between the 
turbulent boundary layer and the shock, and the presence of separated flow may all contrib- 
ute important effects on the loading. The actual physical problem is governed by nonlinear 
partial differential equations with nonlinear boundary conditions, and defies exact solution . 
Hence, we confine ourselves here to an extension of the linear theory only to account for 
the effect of wing thickness insofar as it produces a nonuniform mean flow and, possibly, a 
small locally supersonic region. This should be a significant improvement over the usual 
idealized lifting surface model. 

One method of improving the linearized theory is to account for the nonlinearity in an 
approximate fashion, based on the concept of "local linearization." In this method, the 
effect of thickness is included only to the extent that the flow remains free of shocks. The 
underlying assumption is that the physical state, usually governed by nonlinear equations, is 
adequately described by related linear equations, with all parameters involved having their 
local values taken to be constants in a limited region. This involves replacement of the 
nonlinear equations with linear equations having variable coefficients. 
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These methods were introduced into aerodynamics in an intuitive way for steady flow by 
Spreiter and Alksne (ref. 31), and have now been rigorously validated by means of the 
method of parametric differentiation by Rubbert and Landahl (ref. 32). This approach sug- 
gests that, in the case of unsteady flow, the calculations can be carried out with sufficient 
accuracy using the linearized equations which contain the local values of the steady-flow 
parameters. In his linearized theory, Landahl (ref. 33) cites evidence for the validity of 
applying the concept of local linearization to the case of unsteady flow. 

Moreover, in the case of supersonic flow, it was pointed out by Ashley (ref. 34) that a 
potentially simple way of improving the accuracy of unsteady-flow calculations is to use the 
linearized velocity-potential equation, but with the Mach number of the undisturbed flow 
replaced by the local value which varies spatially due to thickness, mean angle of attack, 
and/or camber-line shape. This work has been extended in two ways. Sankaranarayanan 
and Vijayavittal (ref. 35) applied Ashley's (ref. 34) approach to flow past delta wings and 
found that the general effect of thickness is to reduce the flutter speed; meanwhile, 
Kacprzynski (ref. 36) examined the three-dimensional effects more fully. The above 
investigations are concerned with studying thickness effects in supersonic flow. Except for 
the work of Andrew and Stenton (ref. 37) relatively little information is available about the 
behavior of a finite wing with thickness in the transonic flow regime. 

It should be pointed out that one still earlier method for analyzing the transonic airfoil 
problem also introduces a concept similar to that of local linearization. Using the theo- 
retically exact hodograph transformation, Garrick and Kaplan (refs. 38, 39) developed 
velocity correction formulas which provide a basis for comparison of corresponding compres- 
sible and incompressible flows. In reference 40, Guderley discusses the extension of the 
hodograph method to the problem of three-dimensional transonic flow around finite-span 
wings. However, the advantage of linearization is lost in the extension to three dimensions, 
and his results are limited to certain trends for the influence of aspect ratio for a special 
type of finite-span wing. Guderley concludes that to find a complete solution of the non- 
linear equation for transonic flow around a wing of moderate aspect ratio would prove 
extremely cumbersome, if not impossible. 


Local Linearization Procedure 


The local linearization procedure starts with the unsteady transonic flow equation, 
equation (20), namely 

2 

cp + cp - M (2cp + cp ) = 0 , 
yy zz ^xt Mr 


(42) 


•where the local Mach number M is given by 

1 


M = M 




(43) 
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Here, C pi is the pressure coefficient of the wing associated with the steady flow at its mean 

position. M or C p j will have to be provided either from theoretical results or from available 

test data. The quantity M is all that is required to proceed with the sonic -box method 
solution (ref. 5) . 

For oscillatory motion, equation (42) can be written as (see, e.g. eq. (21)) 

2 

cp + cp - M k(2icp - kcp^) = 0 . (44) 

yy zz x 


Then, by introducing the following modified Prandtl -G lauert transformation 

y = M(x,y) • y 
z = M(x,y) • z 


(45) 


into equation (44) and neglecting the spatial derivatives of M (here M is considered as a 
parameter), one obtains the following 


cp + cp -2ikcp + k cp = 0 
o ~~ - — o o 

y y z z x 


(46) 


where 


cp Q (x,y,z) =M(x,y) cp o (x,y,z). 


(47) 


This equation is in the form that Rodemich and Andrew (ref. 5) used in their sonic-box 
method. Thus a wing properly transformed according to equation (45) may be treated by the 
sonic-box method in the transformed space. Even though the local Mach number may be 
different from unity on a given wing, the consideration of M here is not thought to be the 
primary effect which changes the type of flow; rather, it is considered to be a parameter, 
with the local disturbances unable to travel upstream for the cases under study by virtue of 
the Mach freeze . 


The corresponding linearized tangency condition becomes 


cp - h + ikh 
o~ o o 

z x 


where 


h 0 ( x,y ) = h D ( x 'y) • 


(48) 


(49) 
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Limitations of the Proposed Transformation 


The sonic -box method developed by Rodemich and Andrew (ref. 5) is restricted to a 
wing where no part of it is downstream of any part of its wake(s) . That is, the slope of the 
leading edges of the wing ds/dx is not allowed to change signs as shown in figure 2. The 
leading edge shown in the positive y" quadrant is unacceptable because its slope changes 
sign between the apex and the trailing edge of the wing. While the slope of the leading 
edge shown in the negative y quadrant does not change sign, thus it is acceptable. Since 
the sonic -box method is applied to the transformed wing in the present study, this restriction 
may be expressed as 


i [±7(x) ] * o , 

or, written in the quantities for the physical plane, 

iUsMl 

±s(x) = - M(±s) • .22 , (50) 

where M(±s) represents the local Mach number along the leading edges. 

In addition to the limitation shown in equation (50), there is another limitation con- 
cerning the spanwise Mach number distribution . That is, at x = constant, 

^>0 for 0 < |y| < s . 

a y 

Figure 3 illustrates the nature of the acceptable and the unacceptable functional forms. 
In the sketch, the relation y = y • M(x,y) is drawn as a solid line. As can be seen in the 
illustration, the transformed local span in y ^ 0 is sing I e-va lued (acceptable) whereas that in 
y < 0 is multi-valued (unacceptable) . The multi-valued transformation creates a type of 
fold-over wing in the transformed space which makes this method inapplicable. This condi- 
tion may be relaxed somewhat if at those points where the values of dy/dy = 0 occur within 
the range | y"| < \T\ , y"only deviates slightly from the required monotonica I ly increasing 
nature as y is transformed into y from 0 to s. 
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SONIC-BOX METHOD 


Although any computational technique that is suitable for use with the linearized 
transonic theory could be applied in combination with the present method, the basic com- 
puter program chosen for use in this study was that originally developed by Rodemich and 
Andrew (ref. 5) and known generally as the "sonic -box" program. Since equation (46) and 
the boundary condition, equation (49), to be treated in the transformed space are in the 
form that was used in reference 5, the principal task required to account for the thickness 
effect is the transformation of information between the physical and the transformed spaces 
according to equations (45), (47), and (49). The mathematical development of the numerical 
procedure of carrying out the integration will not be duplicated completely here; however, 
the pertinent equations are included for the purpose of continuity. Since the wing to be 
considered is in the transformed space, the quantities in the transformed space rather than 
that of the physical space are used in the description of the sonic -box method. 


Solution Procedure 


The basic solution satisfying equations (20) and (29) for a thin wing having its mean 
position lying In x^-plane (see fig. 4) representing a point doublet with its axis oriented 
in T-direction is 


0 , x * 0 


= ^ 

1 

ik 'z 


ik | 

f ~2 2\| 

, y + z i 


2rr ~2 ’ 

X 

w 

exp 

" 2 | 

X + - 1 

l X / 


(51) 


, x > 0 . 

Then the solution to equation (20) for a distribution of doublets can be written as 

0 , x 0 


<P 0 (*,y, z ) = 


f f ' ^ x " 5, * )d?dr) ' x > 0 

x-?>0 


(52) 


where p(x,y), representing the doublet strength, may be any function such that the integral 
exists . 

The downwash for a point (xj'y) on z- 0 of the wing mean plane, may be obtained 
from equation (52) 
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yypC?, 1 !) • ^(x-?, Y -“H, z)d?dr| 


w(x,7,0 + ) = Jim ^ W Q \ - J? 

z ~*0 z ~*0 


or 


w 


here 


w(x,7,0 + ) = ff p(? f iD • Kx-?, y-T|, 0 + )d?dri 

<~Sr< 

w+ s 

'Kx,y,0 + ) = lim [cp d (x, y/ z)] 
z'-O 

0 , x £ 0 


= < 


ik ' 1 / y 2 \ 

rST' ex P['2 ,k 


, x > 0 


and 


P(x,y) = cp (x,y,0 + ) 


Thus the boundary value problem becomes 

w(x,y,0 + ) =Jf cp o (?,ri,0 + ) • Kx-5, y-^d^dri, for (x,y) in S , 
and from equation (34) 


(4r* «■) 


cp Q (x,y,0 ) = 0, for (x,y) in W 


( 53 ) 


(54) 


(55) 


Doublet Box Method 


The amplitude of oscillation of the moderately thin wing under consideration is assumed 
to be small. The mean position of the wing is considered to lie approximately in the xy'- 
plane with its nose at the origin and with an unswept trailing edge. The freestream is 
assumed to be parallel to the x-axis and at a Mach number of unity. The value of the 
unsteady potential on the wing may be found by using equation (54). 


To get an approximate solution of equation (54), let the xy-plane be covered with a 
grid of square boxes with the box edges parallel to the x- and y-axis. Let the region E be 
composed of all boxes whose centers lie in S. Thus, E is an approximation to S by boxes. 
The potential cP Q at each box is assumed to be constant, i.e., cp Q . . = constant at(i,j)-th 

box Ejj. The tangency condition is applied to the center (xj, yj) of each box E;j in E and 
the region of integration is replaced by E. Thus equation (54) gives a system of linear 
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algebraic equations for the Sp Q *s: 


M 



• 1*1 

•I 


// 


Y(x. - ?, ^ - Tl)d§dil = w(x. 




(56) 


or 

A(o, |i-j'i)^ 0 =w(x y) - X 2^'-''' IH’I)? ( 57 ) 

°ii' ' 1 i'<i i' i'i' 

where 

A<i — i * , li-j 1 |) = ff Y(x -5, y.-iDd5dn. (58) 

F ' 

i'i' 

If the wing is symmetric about the x-axis ; then only modes of oscillation that are sym- 
metric or antisymmetric in y need be treated. Thus one only needs to consider one-half of 
the wing in a computation. For symmetric modes, equation (57) becomes 



X l>(0, li-i'l) + A(0, i+i'-l)]cp o 

i'*i 1 


= w(x., y.) 


- X Z CA(i-i' , li-j'l) + A(i-i', |j+j'-l|]cp o 
i'<; i'^i 


^ l 


(59) 


The equationsfor j ^0 are implied by j ^1 . 

For antisymmetric modes, equation (59) applies, with the sums of values of A replaced 
by differences . 


DESCRIPTION OF THE COMPUTER PROGRAM 


A major part of the computer program used in this study is adapted from that developed 
by Rodemich and Andrew in reference 5. The capability of the original sonic -box computer 
program is retained and, in addition, it has been extended to include nonlinear thickness 
effects. This is accomplished in an approximate manner by applying the local linearization 
concept in conjunction with the mean local Mach number distribution which is assumed known. 

The wing with thickness is transformed into a wing without thickness in the transformed 
space using the relations in equation (45). The equation to be solved and the boundary 
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condition to be satisfied in the transformed space become those shown in equations (46) and 
(49), respectively. Then the sonic-box method is applied directly, since the linearized 
theory is applicable in the transformed space. 

The leading edge of the planforms of both the physical and transformed wings must not 
have any local maximum in the spanwise direction within the range 0 ^ x < b; i.e. f the 
local semi-span must not decrease as x increases from the nose to the trailing edge. The 
leading edge of any given wing planform is represented by a finite number of straight line 
segments and, based on this, the area is then approximated by a grid of square boxes. In 
each box, represented by its center, the flow parameters are considered to be constants equal 
to their values at the box centers. The trailing edge of the wing must be unswept*, and 
there must be no control surfaces. The computer program as presently formulated considers 
only the symmetric mode* of oscillation, having the xz-plane as the plane of symmetry, 
and the local Mach number at a corresponding point on the upper and lower wing surfaces 
is taken to be the same. 


The doublet potential at the center of each box in the transformed space is obtained by 
using equation (54), with a known down wash, from the sonic-box computer program. Then it 
is converted into the doublet potential of the physical wing by applying equation (47). That 
is, the doublet potential at the center of each box in the physical space is found by locating 
the corresponding point in the transformed space and applying the inverse transformation. 

The inversion results in a non-uniform array, but this method applies interpolation with a 
polynomial fitting procedure to the data before further usage for integrations. Equation (33) 
is used in the calculation of the unsteady pressure coefficient and the generalized aerody- 
namic force coefficients are computed from the following expression (ref, 2): 


L vrwf/K + '^i f j ^ - 


( 60 ) 


where the oscillatory displacement in the z-direction from the mean position of the wing is 
written as 

h 2 (x, y< 0 = (x, y)e lkt = 6.f.(x, y)e‘ kt 


so that the mode shape can be related as f. (x, y) = 


I I 

h Q (x, y) 


(61) 


with h Q as the nondimensional amplitude in z-direction, and 6j as the nondimensional 
amplitude of the j^mode. 

The stability derivatives are calculated by using the following relationships: 


* This restriction is carried over from reference 5 and is not an essential limitation of the 
present method. 
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(62) 


C 


"a 


-ilm CL,,] 


C M- +C M k lm ^ L 22 a ^ L 21 + L 12^ + a L ll^ * ^ 

a q 

The standard output of the computer program consists of the generalized aerodynamic 
force coefficients. In addition to these, the Mach number, downwash, potential, as well as 
the lifting pressure distributions can be printed out if they are desired. 


RESULTS AND DISCUSSION 


Typical unsteady aerodynamic results have been calculated, using the above local 
linearization procedure, for a variety of wing planforms and thicknesses. Such studies are 
necessarily limited by the availability of steady-state data on local pressure or Mach number 
distributions for any desired wing shape or thickness, since this is required as input to the 
basic calculations. The planforms selected for this study are rectangular, delta and parabolic, 
and the steady state Mach number distributions over these wings in a sonic flow field are avail- 
able from reference 41 . This information is furnished to the computer in the form of pressure 
coefficient or Mach number distributions at a group of selected points on the wing. The data 
are then fitted with a pre-selected polynomial form representing the Mach number as a func- 
tion of x and y using the least-square method available in the original sonic-box program. 

The standard output from the analysis is in the form of amplitude and phase of the 
generalized aerodynamic force coefficients as a function of frequency. Typical results are 
shown in figures 5 and 6 for rectangular and delta wings, respectively. The rectangular 
wing has a biconvex airfoil section, with an aspect ratio A =2.0, and is pitching about an 
axis along the leading edge. Based on a distribution of 25 boxes along the centerline chord, 
parts (a) and (b) of figure 5 indicate that for sonic flow conditions, quite large reductions 
in generalized force coefficients are possible for rectangular wings of thickness t = .0521 . 
However, changes as great as 25 percent are larger than anticipated, especially at the 
higher frequencies. In the absence of directly correlated tests, with finite thickness rectang- 
ular wing models, these results should be used with caution since some limitations of the pro- 
gram are more pronounced for this case. The steady Mach number distribution is quite rapidly 
varying, particularly near the leading edge of this wing (ref. 41), so in the transformed 
coordinates the planfbrm is severely distorted. The box size used may be too large to proper- 
ly define these variations. Also, the requirement of the linearized theory that the reduced 
frequency satisfy the relationship k » is probably violated by the large velocity dis- 
turbances near the leading edge. This may explain the apparent lack of convergence in 
magnitude for the two ( L] coefficients in figure 5(a) at frequencies near k = 1,0. At 
high frequency, the linearized theory is expected to be accurate, independent of thickness 
considerations. However, the residual differences between the two results for the higher k 
values used here may exist because that maximum k value is still insufficiently high to insure 
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convergence. This possibility may best be seen in figure 5(b) where good convergence is 
exhibited by the phase angle ©22/ ar) d | 1-22 ! s ^ ows some asymptotic tendency at the highest 
frequency calculated. Further study of the effects of box size is needed for rectangular wing 
cases. 

A typical variation of the generalized aerodynamic force coefficient for pitching moment 
due to pitch oscillations of a delta wing configuration of aspect ratio A = 1.5 are shown in 
figure 6. A distribution of 40 boxes along the centerline chord was used in the computations. 
This wing has an elliptic cross-section in a plane perpendicular to the chordwise axis, and it 
is assumed to be pitching about the apex, or the most forward point of the delta wing; i.e., 
a - 0. Lift and pitching moment magnitudes and phase angles show very little change with 
thickness increase as compared with the rectangular wing values. At medium-to-high reduced 
frequencies, increases in absolute values are of the order of 5 percent, with the maximum 
occurring for 10 percent thickness. Over most of the frequency range, a further increase to 
15 percent thickness results in a reversal of trend, or less increase in magnitude of the force 
coefficients. As will be discussed later, however, variation of pitch axis to more aft loca- 
tions may affect the percent change markedly for a delta wing. 

The magnitude of the unsteady pressure distributions for this delta wing are presented in 
figure 7 for both pitching and plunging motions at reduced frequency 0.1. The pressure 
variation for a change from zero to 10 percent maximum thickness is not very large except 
near the trailing edge. The changes are generally greatest where anticipated based on con- 
sideration of local flow acceleration due to the addition of thickness. 

Comparisons are made in figure 8 with the low-frequency approximation in Landahl's 
linearized theory (ref. 2, chapter 3) for a delta wing without thickness effect. The 
trends of increasing damping-in-pitch (C/y^ + C/y\q), increasing frequency are qualita- 

tively in good agreement for the higher aspect ratio case. For the lower reduced frequencies, 
larger effects on damping-in-pitch are possible due to changes in either thickness or pitch 
axis location, as shown in figures 8 and 9. The differences in absolute magnitudes for the 
zero thickness cases, particularly at higher frequencies, may be rationalized by considering 
limitations of both linearized methods compared in the figures. Landahl’s approximate low 
aspect-ratio method should be adequate at moderately low frequencies for these two cases 
(A = 1.5 and A - 3.0). His solution has an upper frequency limitation, and this may ex- 
plain part of the differences in figure 8. However, the local linearization method also has 
a severe limitation since the potential differences between the "locally constant" function 
over each of the boxes do not cancel at the edges of adjacent boxes when the integration is 
carried out around the boxes. This implies some error in the surface integral which increases 
with the size of the box used. This may explain the lack of convergence with Landahl's 
linearized theory at the higher frequencies. 

As mentioned previously, the delta wing thickness does not appear to affect the results 
as much as for the rectangular wing at the higher reduced frequencies. It is believed, how- 
ever, that some additional effects would be seen if an improved fitting procedure were to be 
employed. The steady flow results used (ref. 41) indicate a less rapidly varying func- 
tion over the surface of the delta wing, but the least-square fitting method in the original 
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sonic-box program essentially reduces the pre-determined polynomial into a linear form. 

The more rapid variations at the leading and trailing edges were ignored. Thus, the trans- 
formed planform was only slightly different from that for the delta wing, and only small 
thickness effects were therefore observed. Further study is planned using improved fitting 
techniques to determine more accurately the sensitivity of the program to delta wing 
thickness variations. 

Only a small amount of*test data are available for oscillatory transonic flow with these 
configurations. The data scatter is generally great, but the more reliable appearing data at 

= 1 in references 42 and 43 are shown in figures 9, 10, and 11. The data for delta 
wings with finite thickness indicate that both trend and magnitudes of the stability derivative 
shown in figure 9 may be in reasonably good agreement at intermediate reduced frequencies. 
The chordwrse pressure distributions plotted in figure 10 show better agreement at the outer 
wing stations. Data from the same tests (ref. 43) in figure 11(b) also show good 
agreement with the analysis for lift-curve slopes as a function of frequency. Trends with 
wing thickness are not observable because of scatter in the test data. 

Only results based on thin-wing theory are available for comparison with the rectangular 
and delta wing calculations in figures 11(a) and 11(b), respectively. Landahl's low- 
frequency approximation (ref. 2 , chapter 4) compares quite well for the rectangular 
wing in figure 11(a) at low reduced frequencies. Similarly, Davies' theory (ref. 3) 
agrees very well with results for the delta wing at aspect ratio, A = 1.5, shown in figure 

11(b). 

The capability of the program to represent other planforms is demonstrated in figure 11(c), 
where a six-segment approximation to the leading edge of a parabolic wing and a distribution 
of 40 boxes along the centerline chord were used. The cross-section consists of a double 
wedge in a plane perpendicular to the chordwise axis. For this case, as for others mentioned, 
the polynomial representing the Mach number distribution on the wing obtained from the least- 
square method was not believed to be of adequate accuracy. Improved procedures are avail- 
able for fitting such functions, and it is recommended that the program be checked for each 
case until modifications can be made in the future. 

Another question concerning the sonic-box procedure in present use is concerned with 
the convergence of the magnitude and phase angle as related to the number of boxes used 
to represent the planform. A typical study which should be made for each new configuration 
is shown in figure 12. It is apparent that the generalized force coefficient converges ade- 
quately for the 40 boxes employed in the delta wing cases studied herein. However, at 
higher frequencies it may be desirable to use a larger number of boxes to properly represent 
the wing planform. For the case with thickness effect, this may be especially needed because 
of the greater complexity of the leading edge of the transformed wing. 

Because of the transformation technique adopted, certain types of wings for which a one- 
to-one transformation cannot be made must be excluded from treatment by the present approach. 
However, this difficulty can be alleviated by using a locally varying source strength distribu- 
tion, related to the local Mach number, in place of the transformation presently employed. 
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(See Addendum to this report.) 


Quantitative evaluation of the local linearization procedure introduced here for unsteady 
flow cannot be completed until other theories with transonic thickness effects are developed, 
or until more reliable test data become available for comparison. 


CONCLUDING REMARKS 


The local linearization concept is applied to the determination of approximate wing 
thickness effects in unsteady transonic flow using the sonic-box computational procedure. 

The principal argument in this concept is that the mean-flow properties in a sufficiently small 
region on a wing may be treated as essentially constant. Thus, the governing potential-flow 
equation becomes linear with variable coefficients. A coordinate transformation then reduces 
this equation to linear form with constant coefficients, which is in the same form as the well- 
known linearized transonic-flow equation. Hence, numerical results in the transformed space 
may be obtained by any computational technique that is suitable for use with the transonic 
linear theory. For demonstration purposes herein, the well-known sonic-box computer pro- 
gram is adopted. 

Sample calculations for delta, rectangular and parabolic wings are made to demonstrate 
the capabilities of the program and to show the contribution of wing thickness effects to 
various unsteady flow quantities. Comparisons are made with results from linearized theories 
as well as with available test data, and reasonable trends with thickness effects are indicated. 

The results for. low aspect-ratio delta wings, having elliptical spanwise cross-sections and 
pitching around the apex, indicate that thickness causes rather small changes in pressure, 
damping-in-pitch, and generalized force coefficients relative to those calculated from linear- 
ized theory using the sonic-box method. However, for this particular wing, these quantities 
do not vary monotonically with thickness; so this trend reversal needs further study with other 
cross-sectional geometries and other pitch axis locations before firm conclusions can be made. 
The zero thickness results for both the delta and the rectangular wings are in reasonable agree- 
ment with the corresponding theory by Landahl for low aspect ratios over the lower and medium 
range of reduced frequencies. 

Although calculations for a low aspect-ratio rectangular wing indicated thickness effects 
as large as 25 percent, caution should be exercised in interpreting these results since certain 
limitations of the computer program and the transformation used, in their adopted forms, are 
most pronounced for this configuration. The wing forms which can be treated with the present 
program are limited to those having a mean local Mach number distribution not very different 
from the freestream value. Those wings which cause a multi-valued transformation (i.e., a 
fold-over wing in the transformed space) must be excluded. In addition, only unswept trail- 
ing edges without control surfaces are allowed by the present computational technique. 
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The fold-over limitation and edge conditions associated with the transformation method 
can be removed by using a modified source strength distribution, related to the local Mach 
number, in place of the transformation presently employed to account for thickness effects. 
Also, improved accuracy may be obtained by replacing the least-square surface-fitting pro- 
cedure used in the original sonic-box computer program with better approximations. 

The limited amount of test data available for delta wings in unsteady transonic flow 
show reasonably good agreement with the analysis for finite thickness cases and for variations 
with reduced frequency. Absolute magnitudes cannot be evaluated very well until more 
reliable test data are obtained. 
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ADDENDUM 


ALTERNATE MODIFICATION OF THE SONIC-BOX METHOD FOR 
UNSTEADY FINITE WINGS WITH THICKNESS EFFECTS 

By S. Y. Ruo 


SUMMARY 


An analytical approach is presented for determining some of the nonlinear characteristics 
of the transonic flow equation for finite thickness wings undergoing oscillatory motion at sonic 
flight speed in an inviscid, shock-free fluid. The thickness effect is introduced into the basic 
solution through the variation of the mean local Mach number on the wing, and the flow 
parameters over the wing are considered to be constants within a small finite region. Thus, 
the nonlinear differential equation can be reduced to a locally linear one with variable 
coefficients. 

The equations required to obtain numerical results are reduced to forms similar to those 
used in the well-known sonic-box method computer program. Thus, with some minor modifi- 
cations to the original computer program, the transonic flow aerodynamic force coefficients 
with thickness effects can easily be obtained. 


INTRODUCTION 


It is well known that a closed-form solution to the small perturbation equation governing 
unsteady transonic potential flow over a finite wing is not easily attainable. The previous 
developments for dealing with the problem of transonic flow over an oscillating wing have 
largely been restricted to linearized theories (see e.g. refs. 2 and 3). These methods 
are a result of uncoupling the unsteady motion from the steady part and are generally limited 
to high frequency cases. Since the steady part is completely uncoupled from the unsteady 
equation, the thickness effects cannot be accounted for by using the linearized theory. 

In the work by Kacprzynski, et al . (ref. 4), it has been indicated that thickness, 
which enters into the mathematical nonlinearity, can significantly affect unsteady pressure 
in the supersonic range. More recently, efforts to approximately incorporate the thickness 
effect in unsteady transonic flow analysis have been made by several investigators (refs. 

28, 29, 30, 45, and 46 ). Among the approaches used are extensions of the concepts lead- 
ing to the parabolic and local linearization methods introduced for steady transonic flow, 
respectively, by Oswatitsch and Keune (ref. 19) and by Spreiter and Alksne (ref. 31) . 
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In the analysis based on the parabolic method concept (refs. 28, 29, 30, 45, and 
46), the thickness effect is considered in a global sense and is determined by the overall 
body shape at its mean position. This approach can be improved by modifying the global 
thickness effect to a locally dependent one to refine the physical interpretation of the thick- 
ness effect in the unsteady transonic flow analysis. In the main part of the present work, 
which is based on the local linearization concept, the thickness effect is considered in terms 
of the mean local Mach number. That is, instead of the freestream Mach number, the mean 
local Mach number is used in the governing equation, which is then a linear equation with 
variable coefficients. With an appropriate coordinate transformation in conjunction with the 
local linearization concept, the equation is converted to one with constant coefficients, and 
the well-known sonic-box method (ref. 5) is used in the transformed space to obtain 
the numerical results. During the course of this study, it was found that the transformation 
used was not always adequate, especially when the Mach number variation over the wing is 
large. In such a case, a one-to-one transformation may become impossible, thus limiting 
the capability of the above technique to account for thickness effects in unsteady transonic 
flow analysis. In both approaches, the parabolic and the coordinate transformation methods, 
the thickness effect is considered to be independent of the body unsteadiness, and only the 
steady conditions at the mean wing position are taken into account. Additionally, the 
solutions are limited to Mod = 1.0, zero mean angle of attack, shock-free attached flow, 
and wings with trailing edges which are unswept and without control surfaces. 

The purpose of this Addendum is to investigate the possibility of formulating the problem 
directly from the doublet potential expressed in terms of the local Mach number and, hence, 
avoid some of the limitations of the transformation method which were discussed in the main 
portion of this report. 

Symbols, references, and figures are as in the main portion of this report. 


EQUATIONS AND SOLUTION 
Governing Equation 

The small perturbation equation governing unsteady transonic flow as considered here for 
M_ = 1.0 in dimensionless form is 

CO 

cp + cp - 2 M 2 cp t - M 2 cp = 0 (A-l) 

T yy ^zz xt T tt 

and for harmonic motion, equation (A-l) becomes 

foyy + <Po„ - 2 *M 2 «> 0x + k 2 M 2 (p 0 = 0, (A-2! 
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where 


_ ikt 

cp - cp Q e . 


Equation (A-2) is a linear equation with variable coefficients. The local Mach number 
on'the wing at its mean position is independent of the motion. 

The coordinate system used is shown in figure 4. 


Boundary Conditions 

The boundary conditions to be satisfied are: (1) any disturbances must vanish at infinity, 
and (2) the flow must be always tangent to the solid surface. The first condition is satisfied 
by the type of solution chosen, and the second condition can be written as 



2 / 0 


0 . 


(A-3) 


The linearized tangency condition for the unsteady part is obtained as 


cp Q = h 0 + ikh Q on z = 0 (A-4) 

z x 

for harmonic motion. 

The I inearized tangency condition shown in equation (A-4) has been completely un- 
coupled from the steady part, and this implies that the mean angle of attack, camber, and 
twist of the wing are small (fig . 1) . 

For = 1.0, as in the linearized case, it is assumed that no disturbances propagate 
upstream. Even though the mean local Mach number on the wing may be greater than, or 
less than, unity, the implication of the procedure is that the locally subsonic or supersonic 
character of the flow is not the predominant influence on the oscillatory potential. 

To z ( x ' Y' 0+ ) = h o x + Ikh o ' for ( x ' y ) In s 

cpo x ( x ' Y> 0 ) + ikc P 0 ( x ' Y> 0+ ) = 0 / (x, y) in w 


cp Q (x, y, 0 + ) - 0 


for (x, y) not in (S + W). 


(A -5) 
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Doublet Solution 


The basic solution, representing a point doublet oriented parallel to the z-axis at the 
origin and satisfying the required condition at infinity fora small finite region on the wing 
where the value of M is considered to be constant, may be written for equation (A-2) as 


cp d (x, y, z;M) = 


0, x^ 0 

ik zM 2 , 1 r 

2rr 3“ - .k C X + 


M 2 (y 2 + z 2 ) 


(A -6) 


]} /■ x>0 , 


X“ z X 

in which M is regarded as a parameter rather than a flow influencing property. 


Equation (A-6) satisfies equation (A-2) only in a small finite region of the wing, so 
the solution may be considered to be of the locally linearized form. More generally, then, 
the solution of equation (A-2) for a distribution of doublets is given by 


9 


^(x , y, z; M) - JJ\ p(§ , Tl)cp d (x-5 , y-'H, z; M) d§dT] , x> 0, (A -7) 
(x - §)> 0 


where p(x, y) is the doublet strength. It has been shown In reference 5 for M = 1.0 that 


+. 


P (x, y) - cp (x, y, 0 ) . 


(A -8) 


From equation (A-7), after taking the partial derivative with respect to z, one obtains 


<Pc 


(x, y, z) = jy d(? , 1)) t(x - 5 , y-T| # z; M)d§dT| , 
(x -§)>0 


(A-9) 


where 


+ Bc Pd 

♦(x,'y, 0 ; M) = 1 .= 

3z l 2 -0* 


0, x ^ 0 

2 


.2 2 


ik M r 1 ., , M" y . i 

~j exp [ - ^ ik (x + _Z_) ], x > 0. 

X 


(A-10) 


In a manner parallel to that for the linearized tangency condition, equation (A-9) may be 
written as 



(A-ll) 
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I 


where w(x,y) = cp Q (x,y,0 ). 

Z 

The doublet strength p can be determined from equation (A-9), or by 
J'J' cp o (x,y,0 + ) Kx-?/ y-'H; M)d§dri= w(x,y), for (x,y) in S 


s+w 


and 


fe- + ik^ cp o (x,y,0 + ) = 0, for (x,y) in 


W 


For a wing with an unswept trailing edge, equation (A-12) becomes 

J'J' ¥> o (x,y,0 + ) *(x-§ f y-i-1; M)d§dri = w(x,y), for (x,y) in S 


(A-12) 


(A-13) 


(A-14) 


Thus, the doublet potential strength cp can be evaluated with a predescribed downwash dis- 
tribution w on the wing planform. 


DOUBLET SONIC -BOX METHOD 


The concept or scheme for obtaining the approximate solution to equation (A-14) 
described in this section is essentially the same as that used in reference 5. The only dif- 
ference is the inclusion of the local Mach number in the integrand of the present case. 


For a wing lying close to the xy-plane, the planform (S) is replaced by a grid of square 
boxes (E) with their centers lying in S. Within each box, represented by its center, the 
flow properties and downwash are considered to be constant. The box having a center at 
(Xj ,yj) is designated by E-. With this replacement, equation (A-14) yields a system of 
linear algebraic equations for the cp 0 ..'s as follows; 



iKx.-S, y.-ri; M)d?dr| = w(x.,y.) 


(A-l 5) 


where M is considered to be a parameter only. 


Defining A as 
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A = // t(x.-§, y.-T|; M)d§d'H 


E 

u ' f • I 

« I 


ik KA 2 

2tt M 


//: 


1 


E <*i - 9 

• fl • J I 

•I 


^ exp 


1 -L 

- 2 lk 


2 2 
M (y. -ri) Z 






and Introducing 

u = (x - 5)/H 
v = M(y - rj)/H 

rn = li - i'l 

n = i - i 1 

t= kH 

into equation (A-16), one obtains 

A= TS M //-7 exp 


■5 it ( u b 


dudv 


dSdri (A-16) 


(A-l 7) 


where integration limits are 


_x-§ 1 1 _ 1 

U 1 H ~ 2 ( | -')- 2 n " 2 


u 2 _ (' “ i') + 2 


= n + 2 


v 2 = M 


MbM< 

i - i 

>-M 

. m ‘i] 

Mb] 


= M 

[-a 


Since the integration in equation (A-l 7) is identical to that formulated by Rodemich 
and Andrew (ref. 5) , the integration method developed in reference Scan be adopted 
directly. The only difference is the inclusion of the local Mach number in equation (A-l 7) 
which reduces identically to the expression used in reference 5 when M 1 everywhere on 
the wing . 
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In order to maintain interchangeability of the sending and the receiving points and to 
minimize the accumulated error of integration around the boxes instead of the wing plan- 
form, the value of M in the expression v is considered to be that of the arithmetic mean 
of receiving and the sending points, i.e., 

M=l[M(x,y) + M(§,iD]. (A- 18) 

The value of M outside the integration in equation (A-17) may be similarly taken as the 
mean value as shown in equation (A-18); but it may be adequate or even more appropriate 
to let M = M(§,rj), i.e., the Mach number of the sending points. This observation is 
based on inspection of equation (A-15) where it is seen that the local velocity potential is 
modified only by the corresponding local Mach number at the center of the i'j'-th box. 

All the equations required for carrying out the calculation of the transonic aerodynamic 
force coefficients with thickness effects have been reduced to the forms similar to those used 
in the original sonic-box method computer program (ref. 5) . Thus it can be quite easily 
modified to obtain numerical results with thickness effects without using the coordinate 
transformation technique described in the main body of the present report. The method 
described herein has avoided the difficulty of the fold-over wing problem encountered 
previously, but the basic restrictions still apply. For example, the Mach number distribu- 
tion on the wing must not be very different from that of the freestream to enable the local 
linearization concept to be applied; and the mean angle of attack, the camber of wing and 
the amplitude of oscillation must also be small so that the linearized tangency condition can 
be used. The mean local Mach number is cqnsidered to be a parameter, only, which cannot 
alter the wave characteristics of the disturbances in unsteady transonic flow. 


CONCLUDING REMARKS 


The approach described applies the local linearization method to obtain an approxima- 
tion to wing thickness effects as represented by the mean local Mach number in the unsteady 
transonic flow analysis. The basic solution of the governing smal I -perturbation velocity- 
potential equation is modified by the mean local Mach number to account for thickness 
effects . 

The previous method for accomplishing this required a coordinate transformation which, 
for certain wings, resulted in a non -one-to-one wing planform transformation. It is shown 
that this new method for introducing the local Mach number into the sonic -box computer 
program does not have this limitation. 
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The equations to be used in the calculation of the aerodynamic force coefficients with 
thickness effects have been reduced to those forms used in the original sonic -box computer 
program; thus, only a minor modification to the original computer program is required to 
obtain numerical results. 
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(a) . Li ft - and phase angle due to plunge. 

Figure 5. - Generalized aerodynamic force coefficients for a rectangular wing of 
aspect ratio 2.0 with a = 0 and 25 boxes chordwise. 
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Figure 7. - Unsteady pressure for a delta wing of aspect ratio 1 .5 with k = 0. 1 . 
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Figure 9.- Comparison of present analysis for pitching -damping derivatives with 
delta-wing test data at various wing thicknesses. 
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Figure 10.- Comparison of chordwise pressure variations with test data at three 
spanwise stations on a delta wing of aspect ratio 1 .442 with 


k = 0.108. 
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NASA-Langley, 1975 - CR“225 










